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We develop a large scale quantum master equation approach to describe dynamical processes of 
practical open quantum systems driven by both coherent and stochastic interactions by including 
more than one thousand true states of the systems, motivated by the development of highly bright 
and fully coherent lasers in the X-ray wavelength regime. The method combines the processes of 
coherent dynamics induced by the X-ray laser and incoherent relaxations due to spontaneous emis¬ 
sions, Auger decays, and electronic collisions. As examples, theoretical investigation of real coherent 
dynamics of inner-shell electrons of a neon gas, irradiated by a high-intensity X-ray laser with a full 
temporal coherence, is carried out with the approach. In contrast to the rate equation treatment, 
we find that coherence can suppress the multiphoton absorptions of a neon gas in the ultra-intense 
X-ray pulse, due to coherence-induced Rabi oscillations and power broadening effects. We study 
the influence of coherence on ionization processes of neon, and directly prove that sequential single¬ 
photon processes for both outer- and inner-shell electrons dominate the ionizations for the recently 
typical experiments with a laser intensity of « 10^® W/cm^. We discuss possible experimental 
implementations such as signatures for coherent evolution of inner-shell electrons via resonance flu¬ 
orescence processes. The approach can also be applied to many different practical open quantum 
systems in atomic, quantum optical, and cold matter systems, which are treated qualitatively by a 
few-level master equation model before. 

PACS numbers: Subject Areas: Atomic and Molecular Physics, Optics, Computational Physics 


I. INTRODUCTION 

Master equation approach is a standard technique for 
open quantum systems [I[ and a successful theory in de¬ 
scriptions of light-matter interactions, such as in con¬ 
densed matter pWsics 0, chemistry and biology Q, 
quantum optics [J, and ultracold gases The mas¬ 
ter equation is quite general and encompasses various 
physical phenomena, as long as these phenomena share 
common physical mechanisms, i.e. the interplay between 
coherence and dissipation. Traditionally, these kinds of 
systems are treated within the framework of few-level 
models [a-Q. With the development of highly bright 
lasers Q, however, the few-level models lose the possi¬ 
bility for describing complex open systems, since energy 
here is deposited in a broad range and relaxes in a vast 
number of decay channels. In this case, new challenges 
appear for real dynamics of these complex systems, and 
a large-scale simulation is inevitable. Here, one open is¬ 
sue, related to coherence effects in dynamical processes 
of complex systems, is still unknown. 

Coherence plays an important role in describing corre¬ 
lation properties of quantum matters and understand¬ 
ing quantum phenomena including lasing Q, Fano 
shape [l^, superconductivity, superfluidity and Bose- 
Einstein condensate [TTI . , and novel phenomena aris¬ 
ing from quantum optics [13 1 and attosecond physics [T3| . 
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To study these coherence-induced quantum features, dra¬ 
matic success has also been achieved in preparing and 
controlling the coherent dominant systems with sup¬ 
pressed dissipations, such as electromagnetically induced 
transparency in quantum optics [l^ and superfluid-Mott 
phase transition in condensed matter physics [l^ . Re¬ 
cently, exploring coherence effects of complex systems 
has been arguably one of the most important topics 
in physics, inspired by generating X-ray laser with ex¬ 
tremely high brightness and fully temporal coherence [Sj , 
such as the Linac Coherent Light Source (LCLS) [1^ . 
where one opened a new era of exploring the interaction 
of high-intensity X rays with complex systems on fem¬ 
tosecond (fs) timescales. Unfortunately, a typical feature 
of X-ray-matter interactions is the rapid decay processes, 
due to the vast relax channels. As a result, it is inevitable 
to investigate the interplay between coherence-induced 
effects and dissipations in X-ray-matter systems. 

Here, the X-ray free electron laser, with an ultrashort 
pulse duration and a high peak brilliance, provides the 
possibility for the study of physical, chemical and bio¬ 
logical properties which are never accessed experimen¬ 
tally before [iSlI . For example, a series of pioneer experi¬ 
ments |19l - [^ focus on understanding of X-ray-matter in¬ 
teractions, including nonlinear photoionizations, hole re¬ 
laxations, Auger electron distributions, and X-ray emis¬ 
sion spectrums in atoms, molecules and solid materials. 
To simulate these experiments, a semiclassical descrip¬ 
tion of radiation-field coupled systems, Einstein’s rate 
equation approach [2§-[3l|, has recently been extended 
to the regime of X-ray-matter interactions for gaining 
fundamental insight into the fast decayed systems, by in- 
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eluding photoexcitation and ionization, electron impact 
excitation and ionization, Auger decay, and their reverse 
processes. In Einstein’s rate equation model, all the ab¬ 
sorption and emission processes are treated with transi¬ 
tion probability and the coherence between different lev¬ 
els is neglected, which are good approximations for either 
incoherent light fields or dominant decay processes of the 
systemic coherence [s^. The agreement that one finds 
between theoretical predictions and experiments [2^ 
at relatively low atomic number densities, verifies the ca¬ 
pability of the model to simulate X-ray-matter interac¬ 
tions. The physical reason is that the free-electron laser 
beam has limited temporal coherence with each pulse 
consisting of a random number of incoherent intensity 
spikes in a fs duration. 

One big challenge, however, which has not been over¬ 
come yet, is the temporal coherence of X-ray free-electron 
laser generated from self-amplified spontaneous emissions 
from electron beams ^3. One goal of the X-ray free- 
electron laser is to enhance temporal coherence of the 
field, and the quest for improved X-ray sources, based on 
the free electron laser, is being tackled with several tech¬ 
niques [ 334 ^ . In fact, the new X-ray pulse with im¬ 
proved temporal coherence [4l| is more suitable to locally 
deposit energy and prepare electronic states, study dy¬ 
namical properties via photon correlation spectroscopy, 
and image biological specimens and long-rang orders in 
liquid and condensed matters [l^. Recently, one gener¬ 
ates the first successful coherent free-electron laser radia¬ 
tion pulses in the soft X-ray regime, based on the seeding 
experiment [d^, which experimentally provides the pos¬ 
sibility for investigating the interplay between coherence 
and dissipation in X-ray-matter systems. Then the cru¬ 
cial issue, related to an X-ray laser with an improved 
temporal coherence, is how to model the ultrafast dy¬ 
namics of the X-ray-matter systems and understand the 
underlying physics. To obtain a more precise description 
of dynamical mechanics, in general, we need a quantum 
mechanical tool for simulations of its time evolution, such 
as time-dependent Schrddinger equation. In contrast to 
the case of low-Z species in dilute gases [43 - l48l| . time- 
dependent Schrddinger equation for complex systems is 
difficult to tackle directly, due to electron-electron corre¬ 
lations, collision processes, and spontaneous and Auger 
decay processes. Approximations for the X-ray-matter 
systems are inevitable, such as master equation approach. 
In the master equation approach, one can include the im¬ 
portant microscopic processes as much as possible and 
investigate the interplay between coherence and dissipa¬ 
tions in the X-ray-matter systems. Instead of a few-level 
simulations one need to simulate real dy¬ 

namics of the complex X-ray-matter systems based on 
large scale simulations, due to the vast decay channels 
induced by the intense X-ray laser. As far as we know, 
this problem is never explored in the framework of quan¬ 
tum master equation approach before, coherent dynamics 
of the complex systems is still unclear, and it is still an 
open issue whether new phenomena arise from coherence 


effects in the ultrafast decayed systems. 

These questions motivate our study in this paper to 
establish a general method for describing coherent dy¬ 
namics of the X-ray-matter systems in the framework of 
master equation approach. In principle, master equa¬ 
tion approach can deal with different kinds of dynam¬ 
ics of the intense X-ray-matter systems by including the 
microscopic processes due to photons, electrons and en¬ 
vironments. In dilute atomic gases for exam¬ 

ple, one should include photoexcitations and ionizations. 
Auger and spontaneous decay processes, and their re¬ 
verse processes. For dilute molecules [^, additional 
processes, such as photon dissociations, should be taken 
into account. For solid-dense warm and hot matters [T7I| . 
electron impact excitations and ionizations play an im¬ 
portant role, in addition to environmental screening ef¬ 
fects [13 . Therefore, we develop a thousand state master 
equation for describing real dynamics of these complex 
systems induced by an X-ray laser, and discuss possi¬ 
ble experimental implementations such as signatures for 
coherent evolution of inner-shell electrons. Our studies 
will provide the basis for understanding the coherence ef¬ 
fects in X-ray absorption mechanisms at a fundamental 
level. Actually, if all the off diagonal elements of the den¬ 
sity matrix are ignored, our method reduces to Einstein’s 
rate equation approach, which is a semiclassical descrip¬ 
tion for X-ray-matter interactions [23 - f^ . The present 
approach can be applied to many different practical open 
quantum systems, which are widely treated qualitatively 
before by a few-level maste r eq uation approach l49ll . such 
as those in atomic physics [50j, quantum optics [5l| and 
cold matter physics [3, . 



Time 


FIG. 1: Sketch of multiphoton absorptions in neon induced by 
ultra-intense X-ray pulses. Excitations and ionizations from 
the Is shell dominate the photoabsorptions, followed by the 
consequent further outer-shell ionizations (a) and Auger decay 
processes (b), ending up with a highly ionized stage. Coherent 
Rabi oscillations interplay with dissipations in the X-ray-atom 
systems, such as Auger and spontaneous decay processes. 

In the recently typical experiment [^ . however, 
atomic gases are dilute and the influence of electron colli¬ 
sions and photon scattering can be neglected [s^ . In this 
paper, we take dilute atomic gases as examples for dis¬ 
cussing coherent dynamics of the rapidly decayed X-ray- 
matter systems (sketch in Fig. [1]), where there are still 
few studies related to real dynamics of complex atoms 
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based on multilevel master equation approach Q. Here, 
comparisons between master equation and rate equation 
approach will be made to investigate the influence of co¬ 
herence on the dynamical mechanics, as related to the 
ongoing experiments with different temporal coherence. 
In parallel, we also perform an approximate calculation 
based on a degenerate atomic master equation approach 
which is computationally more affordable. 

The paper is organized as follows: in section II we 
give a detailed description of the approach. Section III 
covers our results for coherent dynamics of neon induced 
by an X-ray laser, and discusses Auger processes of neon 
based on thousand state master equation approach. We 
summarize with a discussion in Section IV. 


II. THEORETICAL MODEL 
A. Hamiltonian 


We consider a many-body system, such as dilute 


atomic and molecular gases [20|, [21], [26| , and solid-state 
materials [T7| . coupled with incoherent sources such as 
vacuum, irradiated by a high intensity X-ray laser. Inner- 
shell electrons of these atoms and molecules will be ex¬ 
cited, forming a far-off-equilibrium system and typically 
relaxing in a fs timescale via spontaneous, Coulombic 
and Auger decay processes. These processes compete 
with other mechanics, such as processes including coher¬ 
ent photoexcitations and ionizations. Correspondingly, 
the total Hamiltonian of the X-ray-matter systems can 
be written as 


H = Ha+ Hf + Hi + Nine, ( 1 ) 

where the total Hamiltonian is the sum of the Hamilto¬ 
nian Ha of the many-body system in vacuum, the ex¬ 
ternal field Hamiltonian Hp including the coherent and 
incoherent external field, the laser-matter interaction Hj, 
and the incoherent-field-matter interaction Hinc- 
The Hamiltonian 


N 

Ha = efcAfcfc, (2) 

governs the time evolution of the system in vacuum, such 
as dilute atomic and molecular gases, and condensed ma¬ 
terials, with Akk = \k){k\ and N being the total energy 
levels included. Here, \k) and tk denote the eigenstate 
and eigenvalue of the system, respectively. The Hamilto¬ 
nian of the external field is given by 

Hf = Y^ hw^alui + hwfblbi, (3) 

i i 

where h denotes the reduced Planck constant, di (a|) 
denotes the annihilation (creation) operator that corre¬ 
sponds to the i-th mode of the laser with frequency 
and bi and uj^ denote those for the incoherent field. 


With the semi-classical treatment for the laser field, 
the quantization of the light is ignored. Instead the light 
is considered as a electric field, E{t), which interacts with 
the f-th dipole di for the transition between states \k) and 
\k') to give 


i 


where the transition operator Di = \k){k'\, and Rabi fre¬ 
quency hflAt) = ediE{t) at time t with di = (A:'|di|fc). 
We remark here that we have used the dipole approxima¬ 
tion for the coupling in Eq. ®, and can also trace out the 
external field degrees of freedom within the semi-classical 
treatment indicating the Eq. ([3]) can be ignored in the 
simulations. This semiclassical approximation forms the 
basis of most investigations on the many-body systems 
both coherently and incoherently coupled by a strong ex¬ 
ternal field. 

Eor magnetic sublevels, we rewrite the Hamiltonian of 
the matter-field interaction in a more explicit form 

Hi = Y^^iDjj'.e-^‘^^ + il-C.), (5) 

J,J' 

where = eE{t){J\di\J') and 

D,„„ = (-1)^-"^ ( I 


Here a denotes the polarization of the external laser, and 
the dipole operator d^ describes & J ^ J' transition be¬ 
tween states IJ, Mj) and | J', M'j) with Zeeman substruc¬ 
tures. 

Actually, we can remove the explicit time dependence 
from the Hamiltonian, if Rabi frequency is time inde¬ 
pendent and H ei,w, transforming the Hamiltonian 
of the system to an arbitrarily specified rotating frame. 
The interaction Hamiltonian in a frame rotating at the 
laser frequency lo reads 

Hi = ^^(Z4jjv + H.c.). (6) 

j,j> 

The system can also be pumped by an incoherent field, 
such as the black-body radiation field in a hot plasma 
environment, and it reads 

Hinc — ^ ^ ^ • G^bnEi H.C.), (7) 

i n 


where and e„ denote the directions of the dipole mo¬ 
ment di, the polarization of the incoherent field, re¬ 
spectively. Actually, the contributions of the incoherent 
field can directly be included in the master equation ap¬ 
proach in , as discussed in Sec. imi 
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B. Method 


Now, the next question is how to theoretically simulate 
the time evolution of complex systems in the presence of 
the X-ray laser, described by E g. dll). In contrast to low- 
Z atomic and molecular gases 4^4§|. the X-ray-matter 
system is dominated by ultrafast decayed mechanics in 
the experimental timescales, such as spontaneous decay 
processes, where one loses the possibility to keep track of 
the couplings between the system and the environment 
with infinite degrees of freedom. In solid-state materi¬ 
als, electron-ion, electron-electron and ion-ion collisions 
occur rapidly and randomly, where screening and broad¬ 
ening effects due to the solid-density environment should 
be taken into account in dynamical simulations. Hence, 
a statistical description for the X-ray-matter system is 
needed, and here we study the time evolution based on a 
generalized thousand state master equation approach for 
the reduced density matrix of the system, where the de¬ 
grees of freedom of both the environment and X-ray laser 
have been traced out in a perturbative treatment. The 
generalized thousand-level master equation for a complex 
system, coupled to vacuum modes of the electromagnetic 
field and irradiated by an X-ray laser and incoherent ra¬ 
diation field, reads 


^ = -l[iI,p(t)]+£p(t), (8) 

where p = denotes the reduced density ma¬ 

trix operator of the multilevel system, and Cp{t) = 
E,r,/2[2A/3(t)A^ - D\D^m - P(0AA] with r, de¬ 
noting transition rate for the i-th dipole due to the back¬ 
ground radiation pump, spontaneous, Coulombic and 
Auger decays, photodissociations, and collision processes. 
Here, the first term in the right side describes the coher¬ 
ent dynamics of the multilevel system coupled with the 
laser field, and the second term denotes the incoherent 
processes, such as spontaneous and Auger decays, tran¬ 
sitions due to the black-body radiation field and collision 
processes. The contributions of the black-body radiation 
pump are only nontrivial in warm and hot matters, and 
collision processes can be neglected for the dilute atomic 
and molecular gases since it occurs in a fs timescale be¬ 
ing much shorter than the average particle-collision time. 
Note that the broadening contributions of the incoherent- 
field-system interactions in Eq. © and plasma environ¬ 
ments are included to the second term in Eq. ®, while 
the corresponding energy shifts are incorporated in the 
energy levels of the system in vacuum. We remark here 
that the thousand state master equation has unique fea¬ 
tures, which needs large-scale simulations for propagating 
in time a matrix of « 10® x 10® and is not trivial to solve 
directly, such as stability of numerical linear algebra and 
parallel procedure demanding. 

There are two possible ways to take photoionization 
processes into account in our simulations. While we 
treat the ionization as incoherent processes by adding 
photoionization cross section in the incoherent terms for 


the low intensity X-ray laser, we consider ionization as 
coherent processes for multiphoton dominant processes, 
which features analogies to bound-state transitions of the 
system. For multiphoton dominant processes, the ionized 
state composing of the residual system and ionized elec¬ 
trons reads 


1^) — Ijcore, Tj Cfcj ^)i (9) 

where jcorej k, J, e^, P denote the angular momentum 
of the residual system, relativistic angular momentum 
of free electrons, total angular momentum, total energy 
and parity of the system, respectively. In the physical 
systems, these states can be populated by multiphoton 
excitations. Considering the large amount of continuous 
states, selection rules should be used to solely include 
dominant states, such as degenerate initial, intermediate 
and finial continuous states connected by multiphoton 
energies, and at the same time neglect states detuned 
from resonance excitations since the finite time dura¬ 
tion of the X-ray laser implies a finite Fourier width for 
the bound-free or free-free transitions. These selections 
are a good approximation for hydrogen in the presence 
of strong laser beams, and therefore we anticipate that 
the dominant states for ionizations of a complex system 
should be similar. 

Complementary to the multilevel master equation, in 
this work we employ a degenerate master equation ap¬ 
proach to explore the physics of Eq. o, which is compu¬ 
tationally more affordable. Here, the Rabi frequency is 

{2Mj + l){2Mj + l)Oi and the decay rate is defined 
as a total transition probability from one upper state 
I J, Mj) to all the degenerate lower states of the level | J'). 

a. population distribution — Free states \k) form a 
complete basis set, and we can work in the set and obtain 
a number of algebraic equations for Eq. (|5]). Due to infi¬ 
nite number of matrix elements, however, truncation for a 
finite Wstate basis is required. And then the full Eq. o 
is projected onto this subspace spanned by the iV-basis 
states, and it is expected this projection gives the best 
possible description of dynamics of the system induced 
by an X-ray laser, such as a neon gas by including thou¬ 
sands of energy levels. After solving the multilevel mas¬ 
ter equation, we obtain population distributions for each 
level and coherence between different transition states as 
a function of time, by calculating pkk' = {k\p{t)\k'), and 
then comparison can be made with experimental data 
by utilizing ion charge-state spectra recorded by a time- 
of-flight analyzer (20|, even though the recently typical 
free-electron lasers have limited temporal coherence and 
new techniques are needed to improve it. 

Master equation approach can reduce to Einstein’s rate 
equation approach, if the pump field is fully stochastic, 
such as a broadband isotropic light field. In this case, co¬ 
herence embedded in the off-diagonal terms is neglected, 
and the time evolution of the system is characterized by 
the population changes for each energy level. 

b. resonance fluorescence — The spectrum of reso¬ 
nance fluorescence in high-intensity x-ray pulses can also 
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be calculated via master equation approach 
The time-dependent spectrum of the fluorescent light, 
described by autocorrelation function of the electric field 
operator E{t), is given by 

S{t,uj)= [ [ 

Jo Jo 

where E~{t) and E~^{t) denote negative and positive fre¬ 
quency parts of the electric field, respectively. 

The change in time of the occupation number of pho¬ 
tons can be related to the dipole moment operator, and 
it reads 


E-{t) = C{r)Dlit), (11) 

where C'(r) is a proportionality factor at position r and 
can be neglected in a homogeneous system. And then we 
obtain 


/■* 

rt-t2 

/ dt2 

/ Re 

'o 

/o 


■{Dl{h+T)D,{t,)) 




where we only include the contributions in the region 
ti > t2 and introduce the time delay t = ti — t2- 

Thus we need the knowledge of two-time expecta¬ 
tion values of = Dj{ti)Di{t 2 ) for the time- 

dependent fluorescent spectrum of the f-th dipole transi¬ 
tion. Applying the quantum regression theorem , 

we obtain the coupled equation for the i-th dipole tran¬ 
sition, 


— A{ti,t2)] + d^A{ti,t2), (13) 

citi h 

where ti > t 2 , and the initial condition can be ob¬ 
tained by the one-time evolution operator Aij{t 2 ,t 2 ) = 
{Dj{t 2 )Di{t 2 )) via Eq. ([5]). 


C. Examples: dilute atomic gases 

Before proceed to construct the master equation, we 
should first calculate the wavefunctions and energy lev¬ 
els of the system, and obtain the required data, includ¬ 
ing oscillator strength, dipole moment, Rabi frequency, 
spontaneous and Auger decay rates, and photoionization 
cross section. Actually, dilute atomic gases are widely 
used in the recently typical experiment [2^. From now 
on, we take dilute multielectron atomic gases as examples 
for discussing coherent dynamics of the rapidly decayed 
X-ray-matter systems. Without loss of generality, these 
discussions can be easily applied to other X-ray-matter 
systems, such as molecular gases and solid-state mate¬ 
rials by calculating corresponding wavefunctions, energy 
levels and microscopic transition processes. Here, we can 
map the dilute atomic system into a single-atom problem 
each of which can be written as in Eq. ©• The compu¬ 
tations of various atomic radiative processes involve the 


bound and continuum states of different successive ion¬ 
ization stages in the single atom. For the bound states, 
a fully relativistic approach based on the Dirac equa¬ 
tion is utilized, while for the continuum processes the 
distorted wave approximation is employed. The bound 
states of the atomic system are calculated in the config¬ 
uration interaction approximation. The radial orbitals 
for the construction of basis states are derived from a 
modified self-consistent Dirac-Fock-Slater iteration on a 
fictitious mean configuration with fractional occupation 
numbers, representing the average electron cloud of the 
configurations included in the calculation. The detailed 
discussion can be found in Refs. (dOl - f^ . 

c. Rabi frequency — In the dipole approximation, 
the Rabi frequency for the bound-state transition can be 
obtained from the degenerate emission oscillator strength 
of the electric dipole di, 

s/. = 

= 8.749 X 10^® X A£;|(di)|2, (14) 

and then it yields. 

eE - 

= ^{d^) 

= 1.409 X 

^ ( -Mj q M'j^ ^ 

where laser intensity I and transition energy AE are in 
units of W/cm^ and eV, respectively. Here q denotes the 
polarization of the external laser, di denotes the electric 
dipole moment for the J ^ J' transition between states 
\J,Mj) and \J',M'j) with Zeeman structure, rUe denotes 
the mass of atom, E denotes the electric field, and e de¬ 
notes the charge of electron. 

For photoionization processes (bound-free transition), 
the cross section is given by the differential oscillator 
strength 

cr(e) = 

= 8.067 X 10-1® X (16) 

where the energy e is in unit of Ry, a denotes the fine- 
structure constant and oq denotes the Bohr radius. And 
then Rabi oscillation for bound-free processes can be 
written as: 

H = ^(d.) (17) 

= 0.135 xl0i®y g{e)^yi x a{e)/AEde, 

where, I, a and AE are in units of W/cm^, cm^ and eV, 
respectively, and g{€) is the lineshape of the bound-free 
transition, such as broadening due to the finite period of 
time duration for the X-ray laser. 
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For free-free processes, the Rabi frequency of the k —>• 
k' transition with the residual ion in the state |jcore) can 
be given by: 

eE . - . 

^ ^ (jcore; '^Idj I Jcore; ^ ) 5 

h 

where k and J denote the relativistic angular momentum 
of free electrons and the total angular momentum of the 
system, respectively. 

d. Spontaneous decay — Even in the absence of an 
applied field, spontaneous emission cannot be ignored, 
since the excited state interacts with the vacuum fluctu¬ 
ations of the electromagnetic field. After integrating over 
all possible modes and summing over the two orthogonal 
polarizations possible for each wave vector, the sponta¬ 
neous |e) —>■ \g) transition rate is given by a perturbation 
treatment, 


Aj = j g{uj)\n^\'^doj 


=2, ,3 


iireohc^ 


(19) 


where c is the speed of light, eo denotes the permit¬ 
tivity, and cjQ denotes the transition frequency. Here, 
one assumes the lineshape function g(oj) is sharply 
peaked around wq, so that j^g{(jj)\{g\i\e)\^uj^d^uj « 
K5|f|e)Pw^l,5(w)d3w = 27ra;g|(g|f|e)|2. 

e. Stimulated transition — If the |g) —>■ |e) transi¬ 
tion is driven by the external radiation field with photon 
polarization e in the differential solid angle dQ, the tran¬ 
sition rate in this solid angle 0, based on Fermi’s Golden 
rule, is written as: 


J 2heo[2TTcy 

^ith \E\^ = 

As for black-body radiation field, energy density of its 
radiation varies slowly over the range of transition ener¬ 
gies. 


/ g{ijj)p{uj)duj Tsi p(uj) / g{w)dio = 2Trp{uj). ( 21 ) 
J UJ J UJ 


Considering the isotropic and unpolarized properties of 
the radiation field, the total transition rate can be ob¬ 
tained by integrating over all possible modes and sum¬ 
ming over the two polarizations of the field: 


with energy density p{ui) = . We remark here 

that the transition rate Bi = Ain(uj), with n(uj) being the 
mean photon number of the incoherent field at frequency 
UJ [g, lU , such as the black-body radiation field. 

As for laser driven atomic system, one can use the 
relation between laser intensity I and photon number 
with polarization e and wavelength A 

l{Q)/c = nx ,^ 'g fau, (24) 

(ZTTC)'^ 

and then obtains: 

J D r 

id = J 

= J -^^^\{e\r\g)\'^Iie)giuj)duj, (25) 

where one averages over all the polarization directions of 
the dipole moment. Normally, the laser injects from one 
certain direction with central frequency u (I(0)d0 = 
I(u)6(0)d0), and in this case the total transition rate 
can be written as: 

Bge =-^^\{e\r\g)f J I(v)giii)dv (26) 


where g(v) is the line shape of the I 5 ) \e) transition. 

If one defines absorbtion (emission) cross section as: 


ai(v) = Bihv/I(v) 
27r^e^ u 


Seofic 


e|f|5) 5(0> 


one obtains 

H,; = 


Ciiv) ^ dv 


hv 




hv 2 eoTOeC 


f^g(v)dv 


= 109.761 X 10 

f I(u;)f, 


-18 


hAil2i 


huj (huj — hui)"^ + (HAi/2)^ 


(27) 


(28) 


duj. 


where the full-width-half-maximum Av = Ai/27r, and 
laser intensity I and photon energy Huj are in units of 
W/cm^ and eV, respectively. 

/. Autoionization — Complex atoms, irradiated by 
X-ray laser, can form hole atoms and then relax in a fs 
timescale via Auger decay processes, i.e. spontaneously 
emitting one of the electrons in along with refilling the 
hole by outer-shell electrons. The corresponding autoion¬ 
ization rate is expressed as 


and correspondingly Einstein B coefficient 

o2 


^Ein — 


Tre 

3eoh‘^ 


e rl 


( 22 ) 


A? 


2E 



k; j, Mj 


E 

m<.n 


1 


'^mn 



, (29) 


where k is the relativistic angular quantum number of 
the ionized electron, and r^n denotes the relative dis¬ 
placement of between electrons m and n. Here, J and 
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Mj denote total angular moment and magnetic quan¬ 
tum number of the system consisting of the residual ion 
and ionized electrons, respectively. 

For some cases, double Auger processes are nontrivial 
and should be included in the simulations, and the code 
for the atomic data of direct two-electron processes has 
been developed recently [^ . In this work, we do not 
include the contributions from the direct double Auger 
processes, yet its implementation in the master equation 
is straightforward. 

III. RESULTS 

In this section, we discuss coherent dynamics of the 
rapidly decayed X-ray-matter systems. As examples, we 
investigate ultrafast dynamics of inner-shell electrons of 
complex atoms irradiated by an X-ray laser via a thou¬ 
sand state master equation approach. The influence of 
coherence on dynamical mechanics will be investigated 
through comparisons between atomic systems induced by 
X-ray pulses with different temporal coherence. For ion¬ 
izations, multiphoton ejections of electrons from an atom 
subjected to a strong laser field will be investigated in the 
framework of master equation approach, where the reli¬ 
ability of our numerical results will be verified against 
the time-dependent Schrodinger equation. For inner- 
shell dynamical processes, Auger and spontaneous decay 
processes typically occur in a fs timescale and compete 
with other mechanics, such as coherent Rabi oscillations, 
where we will discuss possible experimental implementa¬ 
tions such as signatures for coherent evolution of inner- 
shell electrons. Finally, we will discuss a real coherent 
dynamics of complex atoms induced by an X-ray laser, 
based on master equation approach by including thou¬ 
sands of atomic levels. 


A. Coherent dynamics between bound states 

In this section, we investigate coherent dynamics be¬ 
tween bound states of an atomic system coupled by a 
X-ray laser, and pay special attention to the influence of 
coherence on dynamical evolution by comparising with 
Einstein’s rate equation, which is a reduction of master 
equation approach and a method for the system illumi¬ 
nated by a broadband isotropic light field (incoherent 
light). 

First we study the time evolution of a two-level H- 
like neon with atomic data being obtained via solving 
Dirac equation Due to the pump of the laser 

beam and coupling with the external environment, it is 
expected that the system exhibits a Rabi-oscillating fea¬ 
ture, and then decays to a steady state after a longer 
time. As shown in the upper panel of Fig. [2l we observe 
Rabi oscillations between |J = 0.5, Mj = 0.5) (Is) and 
|J = 1.5, Mj = 1.5) (2p), and, after a few fs, steady 
states are achieved with the excited-state population be- 



t(fs) 



t(fs) 


FIG. 2: Comparison with rate equation. Upper: State pop¬ 
ulations of a two-level H-like neon (between |J = 0.5, Mj — 
0.5) of Is and |J = 1.5, Mj — 1.5) of 2p) as a function of 
time. Lower: Degenerate state populations for the transi¬ 
tion IJ = 0) —^ IJ = 1) (ls^-ls2p) of He-like neon. Here, the 
solid and dashed lines denote populations of the ground and 
excite states, respectively, coupled by a laser with an intensity 
Jo = 10^^ W/cm'^. For the incoherent pump, master equation 
yields identical results with those from the rate equation. 


ing We observe that the stable states 

of master equation at a longer time are identical with 
rate equation approach, whose conclusion is also consis¬ 
tent with analytical results. 

Then we investigate a multilevel atom interacting with 
an isotropic and unpolarized radiation beam, as rele¬ 
vant to Einstein’s rate equation. We focus on the dipole 
transition |J = 0) —>■ |J = 1) (ls^-ls2p) of He-like 
neon, where there are four states in the reduced Hilbert 
space and the upper three levels with angular momentum 
IJ = 1) are triply degenerate due to magnetic splitting. 
We assume the ground state |J = 0) are populated at 
t = 0, and then switch on the pump laser. We observe a 
decayed Rabi-flopping structure in the master equation, 
due to the coupling with the pump laser and the exter¬ 
nal environment, while in the rate equation there is a 
monotonous decay for the ground state, as shown in the 
lower panel of Fig. [51 We also observe that both methods 
cannot yield identical stable states at a longer time. It 















indicates that, for a multilevel atom coupled by a laser 
field, coherence influences both short-time structures and 
long-time stable states, and Einstein’s rate equation ap¬ 
proach cannot give reasonable results since coherence is 
totally neglected in this method. The coherence effects 
can also be clearly seen in Sec. IIII PI 

Another issue needed to be addressed is the influence 
of different temporal coherence of the pump field on the 
dynamics of complex atoms. To simplify our discussion, 
an incoherent thermal radiation field is taken into ac¬ 
count as incoherent pump sources, in addition to a co¬ 
herent X-ray laser. It is expected that the interplay be¬ 
tween coherent and incoherent pump influences dynam¬ 
ical properties of the atomic system. For example, the 
incoherent thermal radiation field can be introduced via 
Rjj' = ryj/h = rj,j-{exp[ft(a;j/ - uij)/kT] - 1}, where 
T is the temperature of the incoherent thermal radiation 
field, n is average number of thermal photons per mode 
at each transition frequency, and Tj j/ denotes the spon¬ 
taneous decay induced by vacuum (the energy shift due 
to coupling with vacuum is included in Hg). In the lower 
panel of Fig. [21 we show results for an atomic system cou¬ 
pled with a coherent field and a thermal radiation field of 
T = 1000 eV, and find that the stable states shift closer 
to those obtained from rate equation. For a fully incoher¬ 
ent pump in the framework of master equation approach, 
we find that the time evolution of the multilevel system 
coincides with Einstein’s rate equation approach. 


B. Coherent dynamics for photoionizations 

g. Comparison with time-dependent Schrodinger 
equation for hydrogen — Ionization is a basic process 
of atoms subjected to an external electromagnetic field. 
For weak laser fields, perturbation theory can be utilized 
to capture these processes. For example, one considers 
a system prepared in an initial state |i) and perturbed 
by a periodic harmonic potential V{t) = which 

is abruptly switched on at time t = 0. Based on time- 
dependent perturbation theory, the first two basis coeffi¬ 
cient c„ can be expanded as: = j^\{f\V\i)\'^S{ujfi — 

uj) corresponding to Fermi’s golden rule, and c„ = 
W Em - 2w) associating with two- 

photon processes, where |*), |to) and |/) denote the ini¬ 
tial, intermediate and final states, respectively. In the 
study of atomic dynamics, one can simply include pho¬ 
ton ionizations as decay coefficients up to the correspond¬ 
ing order, just as used in rate equation approach, where 
the coherence in ionization processes is totally neglected. 
For intense laser fields, however, multiphoton ejections 
of one electron typicall y o ccur, and manifest as domi¬ 
nant mechanics [El. In this case, the perturba¬ 

tion theory lacks the possibility for understanding mul¬ 
tiphoton absorptions, and a nonperturbative treatment 
is inevitable to address this issue. Here, we describe 
coherent multiphoton ionizations in the framework of 
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FIG. 3: Comparison with time-dependent Schrodinger equa¬ 
tion. Ground-state populations of hydrogen as a function 
of time for a laser intensity 7o = 10^“* W /err? with photon 
energies of 14 eV (upper) and 30 eV (lower), respectively, ob¬ 
tained via rate equation, master equation and time-dependent 
Schrodinger equation. 


the atomic master equation approach, and the validity 
of approach will be verified via comparisons with time- 
dependent Schrodinger equation where both the discreet 
and continuous states are treated in the same level. 

In the calculations here, remarkable agreements be¬ 
tween master equation and time dependent Schrodinger 
equation have been obtained for real dynamics of hydro¬ 
gen in the intermediate-wavelength regime. In Fig. [31 
for example, we show the ground-state evolution of hy¬ 
drogen induced by a strong laser field with an intensity 
of 10^"^ W/cm^, and photon energies of 14 eV and 30 
eV, respectively, and find that the master equation ap¬ 
proach can capture photoionization processes of hydro¬ 
gen in the intermediate-wavelength regime. Surprisingly, 
rate equation offers another possibility to address this 
issue at these parameter regimes, even though the de¬ 
tailed structures are different, which indicates that the 
leading order process dominates ionizations here, irre¬ 
spective of coherence or incoherence. We anticipate that 
photoionization mechanism for complex atoms in the X- 
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ray-wavelength regime should also be captured by master 
equation approach. 

h. Coherence in photoionizations of neon — In this 
section, we discuss ionizations of complex atoms in the 
short-wavelength regime. In this regime, time-dependent 
Schrodinger equation encounters difficulties for dealing 
with electron-electron correlations and irreversible pro¬ 
cesses, such as Auger and spontaneous decay processes, 
while rate equation approach, mainly based on photoion¬ 
ization cross sections, fully discards coherence effects in 
the dynamical processes. The underlying physics for ion¬ 
izations of both outer- and inner-shell electrons of com¬ 
plex atoms and its relevant dynamics induced by intense 
X-ray lasers are still unclear. Here, we will address this 
issue related to ionizations of complex atoms driven by 
a strong X-ray laser, based on master equation approach 
whose validity has been verified against ionizations of hy¬ 
drogen as a benchmark. 

In contrast to simple species such as hydrogen, the 
responses of neon subjected to X-ray laser beams for va¬ 
lence and inner-shell electrons are different, since Rabi 
frequency for valence electrons is comparable to free-free 
excitations, while the frequency for inner-shell electrons 
is much larger than free-free processes. We find that 
inner-shell electrons are ionized much faster than outer- 
shell electrons, as shown in Fig. 3] where ionizations of 
valence (upper panel) and inner-shell electrons (lower 
panel) of neon are triggered by an X-ray laser beam with 
a typical experimental intensity of 10^® W/cm^ and pho¬ 
ton energy of 875 eV. Comparisons between rate equa¬ 
tion and master equation are made and remarkable agree¬ 
ments are obtained, even though tiny differences can be 
found in the detailed structures, due to coherent oscil¬ 
lations for ionizations in master equation approach. As 
far as we know, this is the first direct proof for coher¬ 
ent ionizations of both outer- and inner-shell electrons 
of neon in the short-wavelength regime, where sequen¬ 
tial single-photon processes dominate the absorptions for 
the recently typical experimental conditions. This con¬ 
clusion is consistent with the recent X-ray experiments, 
where one did not observe clear signatures of multipho¬ 
ton single-electron ionizations. Hence it is not surprising 
that comparisons between theories and experiments sug¬ 
gest the dominant absorptions of electrons are sequential 
single-photon processes in Young’s experiment [20j. 


C. Coherent dynamics in Auger decay processes 

In addition to ionizations of inner-shell electrons. 
Auger decay is another dominant process for complex 
atoms, which carries valuable information about inner- 
shell electronic structures and dynamical properties of 
atoms, molecules, and solids. Specifically, inner-shell 
electrons can be excited or ionized by the X-ray laser 
beam and a hole atom is formed. The atom is normally 
unstable and relaxes through Auger and spontaneous de¬ 
cay processes, which typically occurs in the fs timescale 



FIG. 4: Ionization of valence electrons (upper) and inner-shell 
electrons (lower) of a neon gas induced by an X-ray laser for a 
laser intensity of 10^® W/cm^ with photon energy of 875 eV, 
obtained by master equation (red) and Einstein’s rate equa¬ 
tion approach (blue), respectively. We find that single-photon 
processes for both outer- and inner-shell electrons dominate 
ionizations for the recently typical experiments with a laser 
intensity of « 10^® W/cm^. 


and can smear out coherence-induced signals embedded 
in Rabi oscillations. In this section, we will investigate 
the competitions between Rabi oscillations, photoioniza¬ 
tions, Auger and spontaneous decays, since the magni¬ 
tudes of these processes are typically of the same or¬ 
der. These processes, however, are normally beyond the 
time revolution of the X-ray experiments. Fortunately, 
resonant fluorescence provides another available tool for 
studying these competitions and underlying coherence ef¬ 
fects. 

In our simulations, we study dynamical evolution of 
ground-state neon subjected to X-ray pulses with dif¬ 
ferent temporal coherence and resonant photon energy 
relative to ls^2s^2p® —>■ ls^2s^2p®3p^ transition. After 
shining the X-ray laser on neon, inner-shell electrons are 
coherently excited and coupled between Is and 3p or¬ 
bitals. On the other hand, the hole state ls^2s^2p®3p^ 
is unstable and decays in a fs timescale, mainly due to 
Auger decay processes by refilling the inner Is orbital 
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FIG. 5: Competition between Rabi oscillations, Auger and spontaneous decays, and photoionizations of a neon gas induced 
by an X-ray laser beam with a peak intensity of 10^® W/cm^. (a) Populations of degenerate states (solid) and 

ls^2s^2p® (dashed) coupled by a flat-topped X-ray laser beam with a resonant photon energy for the Is —>■ 3p transition, 
obtained via master equation (red) and rate equation approach (blue), with the inset for the fluorescence spectrum of the 
ls^2s^2p® —>■ ls^2s^2p®3p^ transition, (b) Ground-state populations of Ne for a resonant Is —>■ 3p (red) and a red detuned 
pump A = —2 eV (blue), where the inset shows the fluorescence spectrum of the ls^2s^2p® —>■ ls^2s^2p®3p^ transition with 
the black dashed line being the power-broadening guideline for the central peak, (c) Neon gases driven by a Gaussian X-ray 
pulse of a FWHM duration 15 fs, with the inset for fluorescence spectrum of ls^2s^2p® —>■ ls^2s^2p®3p^ transition (red) and 
those driven by a flat-topped pulse with the identical fluence (blue). 


by the outer-shell electrons, or spontaneous decays of 2p 
and 3p electrons. The hole state can also be further ion¬ 
ized with double-core forming or via sequential valence 
electrons. To investigate the interplay between these dif¬ 
ferent processes driven by X-ray pulses, we include all 
the dominant microscopic processes. We observe a de¬ 
cayed Rabi-flopping structure between Is and 3p states 
coupled by an X-ray laser, based on master equation ap¬ 
proach, while it is absent in rate equation method with 
a monotonous population changing. It indicates that the 
atomic coherence is distinctly embodied in Rabi oscilla¬ 
tions driven by the X-ray laser, in spite of extremely fast 
Auger decays in the fs timescale. At present, however, 
one lacks reliable tools for observing these fast Rabi os¬ 
cillations by real-time X-ray images experimentally. For¬ 
tunately, the coherent dynamical information for inner- 
shell processes can be extracted from fluorescence spec¬ 
tra, based on master equation approach. 

As shown in the inset of Fig. [Ha), we observe a 
triple-peak structure for the resonance fluorescence of 
the ls^2s^2p® ^ ls^2s^2p®3p^ transition, while only one 
central peak occurs for incoherent light pump. Con¬ 
sidering the remarkable shift (^ 2 eV) of the satel¬ 
lite line from the central peak for a laser intensity of 
10^®W/cm^, the triple-peak spectra provide a valuable 
tool for studying Rabi floppings and are expected to be 
detected via recording photons from spontaneous radia¬ 
tive decays [l^. The next issue is related to the sta¬ 
bility of the triple-peak structure against different ex¬ 
perimental conditions, such as laser detuning and pulse 
shapes. We find that the triple-peak structure is stable. 
Specifically, we observe that the side peaks demonstrate 
asymmetric along with a red shift of the central peak for 
a red-detuned pump, as shown in Fig. IHb), whereas a 
Gaussian X-ray pulse with a FWHM duration of 15 fs 
yeids a broadened triple-peak structure compared to the 
flat-topped pulse, as shown in Fig. IHc). Note that we 


also observe a ~ 5 eV broadening due to the extreme 
strong X-ray laser with an intensity of 10^® W/cm^, as 
shown in the guideline in the inset of Fig. IHb), which 
is normally referred to as power broadening as will be 
discussed in Sec. mmi We expect that our discussions 
provide valuable insight for investigating inner-shell co¬ 
herent dynamics in the upcoming experiments. 

D. Real coherent dynamics of a neon gas irradiated 
by an X-ray laser 
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FIG. 6: Neon charge-state yields by a far off-resonant beam 
with photon energy of 800 eV. Good agreements between dif¬ 
ferent theories [13, [^ and experiments indicate that co¬ 
herence, induced by a far off-resonant laser beam, plays a 
tiny role in the time evolution for the present experimental 
conditions. 

In former sections, coherence effects have been investi¬ 
gated through ultrafast dynamics of model systems cou¬ 
pled by X-ray pulses with different temporal coherence. 


































11 


In this section, we will investigate coherence effects be¬ 
yond model systems, and take neon as examples for real 
dynamics of complex atoms induced by an X-ray laser, 
based on master equation approach by including thou¬ 
sands of atomic levels. The X-ray laser can sequentially 
excite inner-shell electrons and creates a nonequilibrium 
state decayed through a vast number of channels. To sim¬ 
plify our simulations, we add the atomic orbitals gradu¬ 
ally, and mainly focus on dynamics of the lowest-lying 
orbitals of neon, such as Is, 2s, 2p, 3n and 4n, and 
inner-shell excited states ls2n^, Is2n3p and Is2n4p as 
well, which normally includes atomic levels up to an or¬ 
der of 10^ and spans the Hilbert subspace as basis states 
for real dynamics of neon. Our selections are supported 
by the recent R-matrix calculations for inner-shell elec¬ 
trons [^. After obtaining atomic data, we study the 
interplay between coherence-induced effects and dissipa¬ 
tions in X-ray-atom systems, based on large scale simu¬ 
lations. We will verify to which extent coherence demon¬ 
strates a dominant role in the ultrafast decayed nonequi¬ 
librium system. 

i. Off-resonant dynamics of neon in Young’s 
experiment — First we review the details of the X-ray 
free-electron-laser experiment for neon in Ref. (2^. In 
this experiment, X-ray pulses with photon energies of 
800, 1050 and 2000 eV are injected into a neon gas in the 
atomic chamber, and they observe rapid photoabsorp¬ 
tions of atomic gases in the fs timescale. As pointed out 
in the experiment, effects due to particle collisions and 
photon scattering can be neglected, since the neon gas is 
dilute and induced by pulses with photon energies below 
8 keV. Therefore, photoabsorptions are the dominant 
processes of neon triggered by the X-ray beam in the 
ultra-intense, short-wavelength regime. Here, we take 
photon energy of 800 eV as examples for investigating 
the influence of coherence on the sequential multiphoton 
ionizations of the neon gas, where the photon energy is 
far below resonant Is —>■ 3p excitations and it is referred 
to as off-resonant dynamics of neon. Remarkabl e ag ree- 
ments are obtained between different theories |20l . l30| 
and experiments [l^, as shown in Fig.[6l which indicates 
that coherence plays a tiny role in the time evolution 
of neon induced by a far off-resonant X-ray pulse in the 
present experiments. For comparisons, here we use both 
Gaussian and flat-topped pulses to simulate the exper¬ 
imentally accessible dynamics of inner-shell electrons 
of neon, and find that the charge-state distributions is 
insensitive to the pulse shapes. 

j. Near-resonant dynamics of inner-shell electrons of 
neon — In this paragraph, we study the time evolu¬ 
tion of neon subjected to a strong X-ray laser field 
with near-resonant photon energies, based on large scale 
simulations by including energy levels in an order of 
lOA In Fig. [3 we demonstrate charge state popula¬ 
tions as a function of time for a laser intensity Iq = 
2.5 X lO^^W/cm^, which is a typical free-electron-laser 
intensity in the recent experiments, obtained via mas¬ 
ter equation [solid line in Fig. [Zl[a),(b),(c)] and rate 
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FIG. 7: Coherence-induced suppression on time evolution: 
charge-state populations of neon induced by the X-ray laser 
with an intensity of 2.5 x lO^^W/cm^ for different photon 
energies: near-resonant excitation with a red detuning of 15 
eV (a) and resonant photoexcitation [(b),(c)] for ls^2s^2p® —^ 
ls^2s^2p®3p^ transition, and for different pulse shapes: flat- 
topped [(a),(b)] and Gaussian pulses of a FWHM duration 25 
fs (c). Discrepancies between master equation and rate equa¬ 
tion are the results of coherence-induced Rabi oscillations and 
power-broadening effects. Inset: Fraction yields for different 
charge stages of neon obtained by master equation (red) and 
rate equation (green) for a pulse duration of 200 fs (a), and 
induced by a Gaussian pulse (red) of a FWHM duration 25 
fs and a fiat-topped pulse (green) with the same fluence (c). 
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equation approach [dashed line in Fig. [7l[a)]. The pho¬ 
ton energies are chosen near resonant frequency relative 
to the ls^ 2 s^ 2 p® —?> ls^2s^p®3p^ transition, such as a 
red shift of 15 eV [Fig. |71[a),(b)] and the resonant case 
[Fig.|ZIb),(c)], which are both in the photon-energy range 
of LCLS (800-2000 eV). We find that coherence can sup¬ 
press the multiphoton ionizations of neon induced by 
the ultra-intense X-ray pulse. As shown in the inest of 
Fig. Era), for example, the charge state population of 
Ne^^ after a 200 fs evolution, is twice bigger than that ob¬ 
tained by rate equation approach for incoherent pumping 
with an intensity of 2.5 x 10^^ W/cm^. The physical ori¬ 
gin of the discrepancy in the time evolution is the results 
of two effects, which are neglected in the Einstein’s rate 
equation. The first one is the coherence in the inner-shell 
resonant absorption processes, i.e Is —>■ 3p excitations for 
Ne and Is —>■ 2p for Ne^+, Ne'^+ and Ne^’*'. In contrast 
to monotonous changes in rate equation approach, real 
physical processes between different energy levels demon¬ 
strate a Rabi-flopping structure, due to the light-atom 
coupling for the Is —^ 2p and Is 3p transitions. The 
second one is the power broadening effects, due to the 
extremely strong X-ray laser field, which can be up to an 
order of a few eV at a laser intensity of 2.5 x 10^^ W/cm^, 
based on a two-level estimation \/^‘^ + 2^^ with 7 and 

being spontaneous decay rate and Rabi frequency, re¬ 
spectively. This point can also be verified via fluorescence 
spectra, as shown in the inset of Fig. |Sl[b). Interestingly, 
we find that the red-shift case (relative to Is —>■ 3p tran¬ 
sition) ionizes more electrons after 150 fs evolution, since 
the Is —)• 2 p excitations of inner-shell electrons come into 
play for charge states Ne^^, Ne^^ and Ne®^, whereas the 
resonant case ionizes electrons fast in the early stage af¬ 
ter subjected to the X-ray laser field, since the resonant 
Is —3p excitations dominate photon absorptions of Ne 
in the first 50 fs, as shown in Fig. Etb). The influence of 
temporal pulse shape on the charge-state distributions is 
also discussed in Fig. EKc), where comparison has been 
made between flat-topped (green) and Gaussian pulses 
(red) in the inset. We remark here that the temporal 
pulse shape has limited impact on the charge state dis¬ 
tributions, whose conclusion is consistent with those for 
incoherent pulses [ 2 ^. 

Next, we address the issue related to power broaden¬ 
ing effects in more details. Actually, the broadening ef¬ 
fect, resulting in the line shape of a dipole transition in 
an atom, is a basic feature for describing electron mo¬ 
tions and interactions with external fields. The under¬ 
lying mechanisms of broadening effects are diverse, in¬ 
cluding natural broadening due to spontaneous decays of 
excited states, Doppler broadening due to thermal mo¬ 
tions of atoms, collisional broadening due to collisions 
with other atoms or ions, and Stark broadening due to 
energy shifts induced by an external field. Here, we find 
that the dominant broadening effect for the dilute neon 
gas in the LCLS experiment is the power broadening (« 1 
eV) due to the extremely strong external laser, which is 
up to 10 times bigger than those from Auger and spon- 
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FIG. 8: Power broadening effects: Ne fraction in a neon gas 
subjected to a Gaussian X-ray pulse with an intensity of 2.5 x 
lO^^W/cm’^ and a FWHM duration of 10 fs, for different 
photon energies, obtained by master equation (red solid) and 
rate equation approach (blue dashed). 

taneous decay processes (~ 0.1 eV). Normally the rate 
equation approach does not include the power broaden¬ 
ing effects in the present calculations, while in the master 
equation approach they are included automatically. In 
Fig.l we show the Ne fractions subjected to a Gaus¬ 
sian X-ray laser with an intensity of 2.5 x lO^^W/cm^ 
and a FWHM duration of 10 fs, where two local minima, 
obtained from rate equation (blue dashed), correspond 
to resonant excitations for ls^ 2 s^ 2 p® —>• ls^2s^2p^3p^ 
and ls^ 2 s^ 2 p® —>■ ls^2s^2p®4p^ transitions, respectively. 
However, the oscillating structures for resonant excita¬ 
tions at different photon energies are smoothen in the 
results obtained from master equation approach (red cir¬ 
cle), due to power broadening effects. Note here that the 
oscillating structures can occur at other photon energies 
associated with Is —1 np (n > 5) transitions, but are be¬ 
yond the spectrum resolution in the present experiments 
and an improved sharp-band X-ray laser is needed in the 
future. 


IV. SUMMARY AND OUTLOOK 

In conclusion, we establish a general method for de¬ 
scribing coherent dynamics of X-ray-matter systems in 
the framework of master equation approach. We take 
dilute atomic gases as examples for discussing real coher¬ 
ent dynamics of the rapidly decayed X-ray-matter sys¬ 
tem, based on a thousand state atomic master equation 
approach, by including coherent pump and incoherent re¬ 
laxations due to spontaneous and Auger decays. We find 
that coherence can suppress the sequential single-photon 
ionizations of a neon gas in the ultra-intense X-ray field, 
compared to the rate equation approach, and the phys¬ 
ical reason can attribute to the coherence-induced Rabi 
oscillations and power broadening effects, which are both 
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neglected in the Einstein’s rate equation. We also find 
that single-photon ionizations for both outer- and inner- 
shell electrons dominate the absorptions of a neon gas 
for the recently typical experiments with a laser inten¬ 
sity of ~ 10^® W/cm^, irrespective of coherence. A typi¬ 
cal feature of coherent evolution of inner-shell electrons is 
Rabi oscillations with a frequency in the order of 10^® Hz, 
which is beyond the current experimental resolutions in 
the time domain. Instead, we discuss resonance fluores¬ 
cence spectra for possible experimental implementations 
for coherent dynamics of inner-shell electrons. 

With the quick development of free-electron lasers, 
temporal coherence of the X-ray laser is improved ex¬ 
perimentally, which provides the possibility for studying 
time-resolved coherent phenomena in atoms, molecules 
and solid-state materials [b^. There are still a number 
of open issues referred to coherent dynamics of complex 
systems irradiated by an intense X-ray laser. For the 
methods described here, we are able to discuss hollow- 
atom signatures and its corresponding coherent dynam¬ 
ics, where the competition between the double- and 
single-hole generations in coherent evolution is still un¬ 


clear. Moreover, we can also investigate multiphoton pro¬ 
cesses of atoms and molecules in gases and solid-dense 
matters, such as two-photon KK-shell transitions, where 
the issues are related to non-sequential double ioniza¬ 
tions and electron-electron correlations. Finally, in a 
dense hot plasma environment generated by dense gases 
or solid-state materials induced by an X-ray laser, the en¬ 
vironment is much more complicated, both radiative and 
particle-collision processes should be taken into account 
for understanding the phenomena of the system. 
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